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Abstract 

Recent experiments on quantum turbulence in superfiuid helium made use of small tracer par- 
ticles to track the motion of quantized vortices, determine velocity statistics and visualize vortex 
reconnections. A problem with this visualization technique is that it may change the turbulent 
flow which it visualizes. To address this problem, we derive and solve the equations of motion of a 
quantized vortex ring which contains a given number of particles trapped on the vortex core, hence 
derive a quantitative criterion to determine in which measure small particles trapped in quantized 
vortices affect vortex motion. Finally we apply the criterion to a recent experiment. 

PACS numbers: 67.25.dk, 47.37.+q, 47.32.cf 
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I. MOTIVATION 



The study of quantized vortices ^ in superfluid helium (He II) and of quantum turbulence 
, |8j has been held back over the years by the difficulty of flow visualization 
near absolute zero. Fortunately the problem has been recognized, and new visualization 
techniques have become recently available. One of the most promising technique is based on 
trapping small micron-size tracer particles (made of glass, polymers or solid hydrogen) onto 
quantized vortices flfiol . The use of tracer particles has made possible, for example, to study 
velocity statistics [ll| in superfluid turbulence and to visualize 12| individual reconnections 
of quantized vortices, a process which is crucial in the dynamics of turbulence . 

To interpret these recent experiments it is necessary to understand the interaction of the 



tracer particles with the quantized vortices 
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161]. In particular, we need to find in 



which measure the tracer particles disturb the vortices, hence change the flow which one 
tries to visualize. Clearly, on one hand a large number of tracers improves the images and 
the signal to analyze, but on the other hand too many tracers must affect the motion of 
vortices in a significant way. What is the maximum density of tracers which can be used? 

To answer this practical question, we introduce a quantitive measure of the tracers' 
disturbance, which we call the superfluid Stokes number. After deriving and solving the 
governing equations, we determine by what amount the motion of a quantized vortex ring 
of radius R is disturbed by the presence of tracer particles of radius a trapped into the 
vortex core. By interpreting R as the typical local radius of curvature i of the filaments 
(i ^ where L is the vortex line density), we apply the result to a recent experiment 

in which quantum turbulence was visualized. 



II. EQUATIONS OF MOTION 

Consider a quantized vortex ring of radius R which moves in the z direction in superfluid 
helium at temperature T. small buoyant tracer particles of mass m and radius a are 
trapped in the vortex core. For the sake of simplicity we assume that the vortex ring 
remains axisymmetric during the evolution and that the particles remain trapped in it. Let 
(r, 6, z) be cylindrical coordinates. The vortex position and velocity are = (i?, 0, z) and 
^^L = (-R, 0, i) = (m, 0,f), where R = R{t), z = z{t), t is time, and a dot denotes a time 
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derivative. Let the circulation vector he k = (0, k, 0) where k = 9.97 x 10"'^ cm^/s, and let 
k = K,/ K, he the unit vector tangent to the ring in the 6 direction. The forces acting on the 
vortex ring are the Magnus force Fm, the friction force with the normal fluid, and the 
Stokes force F5 induced by the tracer particles. The Magnus force is Fm = 27c RiM, where 



(1) 



Ps is the superfluid density, v^t = + Vj is the total superfluid velocity at the vortex line, 
Vj = (0, 0, Vi) is the self-induced velocity of the ring at T = 0, where 



Vi 



4nR 



In {8R/0 - 2 



(2) 



^ 10"^ cm is the vortex core radius, and = (0, 0, Vs) is an externally applied superfiow, 
which for the sake of simplicity we assume in the z direction. The friction force Fd arises 
from the interaction of phonons and rotons (which make up the normal fluid) with the vortex 
core; we write F/j = 271 Rf£,, where the friction per unit length is 



ffl = Psi^[Toi^n - vl) + T'^k X (v„ - Vl)], 



(3) 



and v„ = (0, 0, f „) is an externally applied normal fluid velocity, which again we assume to 
be in the z direction. The quantities Fq and Fq are dimensionless temperature-dependent 
friction coefficients related to the friction coefficients 70 and 79 calculated by Barenghi, 
Donnelly & Vinen 17|] by 



7o 



F' 



7o 



(4) 



Note that Fq and Fq can be expressed in terms of the more used friction coefficients B 
and -B' as [17 1 

n n R 

(5) 



PnP 

7o = — 



2p [l-B'pJ{2p)f + B^pl/{Ap^y 

, PnP [B^Pn/{2p)-B'{l-B'pr,/{2p)] 

7o = —1^- 



2p'^ [l-B'pJ{2p)f + B^pl/{Ap^) ■ 
If spherical particles of mass m = 47ra'^p/3 each, where pp is the particle density, are 
attached to the vortex ring, the ring has total mass mN and experiences the Stokes force 



Fs = 6TTaUnPnN{Vn " V^), 



(7) 



where Un = p/pn is the normal fluid kinematic viscosity, p is the viscosity, and p„ is the 
normal fluid density. This simple linear Stokes drag is appropriate because the Reynolds 



number based on the particle size is extremely small. For the sake of simplicity and to make 



contact with some experiments 
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12| we assume that the particles are buoyant, i.e. they 



have density Pp = p = Pn + Ps, hence there is no Archimedes force. 
The equations of motion of the vortex ring are thus 



dt 



(8) 



mN-^ = 2TiR{hi + f^) + F5. (9) 

The resulting equations for -R(t), z{t), u{t) = R and v(t) = z are: 

R = u, (10) 

z = v, (11) 

it = cR[v -Vs~Vi- Vqu + r[,(f„ - v)\ M, (12) 

r 

V = cR[-{l - T'q)u + To{Vn - v)] + -{vn - v), (13) 

r 



where 



2ttPsK m 



(14) 



mN GnaUnPn 
The quantity r is the Stokes relaxation time of the particles. 

It is convenient to rewrite the equations of motion in dimensionless form. We choose the 
initial radius Rq = -R(O) as the unit of length, and Vq = k/Rq as the unit of speed; the unit 
of time is thus RI/k. We obtain 

R = u, (15) 

z = V, (16) 

V -Vs-Vi-rou + r'^lvn -v) = ^{eii + (u), (17) 

K 

- (1 - r',)u + Toivn -v) = ^{ev + (v), (18) 

where t, R, z, u, v, Vt, Vn and Vg are now dimensionless. The dimensionless quantities e and 
C are defined as 

2N / 



, , p ; V r C = 3iV p^ - p . (19) 

3 \RoJ \psj \aJ\psJ\KJ 

The above dimensionless form of the equations of motion allows us to compare the motion 
of a vortex ring with tracers against the motion of a bare vortex ring [N = 0). If = 



then e = C = 0, and Eqs. (HID-dlSD reduce to 



V -Vs-Vi- TqU + T'Q{Vn -V) =0, 



(20) 



il-T'^)u + To{vn-v) = 0. 



(21) 




(22) 



[(1 - r'o)^ + rg] 



(23) 



as found in Ref. Il7|. If we further set Vn = Vg = we obtain 



R = u 



(24) 



Z = V = 



(1 - 



which means that, as expected, the ring shrinks {u < 0) and speeds up, transferring its 
kinetic energy to the stationary background normal fluid. If we further set T = we obtain 
the expected stationary solution u = 0, v = Vi. 

III. RESULTS 

We choose values of parameters which are relevant to typical experiments (T = 2 K, 
a = 10~^ cm and Rq = 10^^ cm) and solve Eqs. f|T5|) - f|T8|) numerically using the Adams- 
Bashforth method. In all calculations we take Vn = Vg = for simplicity. The initial 
conditions are given by Eqs. fl22l) and fl23|) . During the evolution the vortex ring shrinks. 
Clearly, when a significant fraction of the circumference of the ring is covered by tracers, 
our model breaks down, as there is not enough vortex length to provide the correct induced 
velocity fj. Since we are not interested in this last stage of the evolution, we arbitrarily stop 
our calculations when the tracer particles cover one third of the circumference, that is, when 
the dimensionless ring's radius becomes smaller than {3N/iT){a/Ro). 

Fig. [T] compares the evolution of the dimensionless radius R{t) of vortex rings with = 
10, 20, 40 and 80 trapped tracer particles to the decay of a bare vortex ring (A^ = 0, solid 
curve labeled "a"), computed solving Eqs. (l22l) and ( l23l) . It is apparent that the more 
particles are trapped (hence the larger e and ( are), the faster the decay of the ring is 
compared to the decay of a bare ring. 
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Fig. [2] plots the corresponding time dependence of the radial velocity component u = R. 
Initially the ring shrinks at constant rate, then the radius decreases faster and faster. 

The axial velocity component i = f is shown in Fig. [3l Initially v remains similar to the 
velocity of the bare ring; during the final part of the evolution a vortex ring with tracers 
moves much faster than a bare vortex ring, because, as shown in Fig. [H the radius has 
become much smaller. 

The results described above suggest that the evolution of vortex rings which contain 
tracers is sufficiently similar to the evolution of a bare vortex ring, provided that the ring 
is not too small and that it does not contain too many tracers. The important question 
is how to quantify the perturbation which the tracers induce on the the ring's evolution so 
that one can decide whether the visualization with tracers affects or not the evolution of 
vortices. To answer this question we start by remarking that the ring's decay is caused by 
both friction and Stokes forces; going back to dimensional variables, we note that the ratio 
of the magnitudes of these forces is 

We call S the superfluid Stokes number. In the absence of external forcing (t>„ = = 0) 
the superfluid Stokes number decreases during the evolution, starting from the initial value 
S'(O) = C/ro at t = 0, because R decreases with time. Fig. H] shows the time evolution of S 
corresponding to Figs. [1], [2], andO It is apparent that if 5* <C 1 the evolution of a vortex ring 
with tracers is very similar to that of a bare vortex ring, but if S becomes of order unity or 
larger then the tracers affect the ring's motion significantly. 

To put in evidence the temperature dependence of the superfluid Stokes number, we 
rewrite Eq. (l25i) as 

S = 6(3, (26) 

where 

- (27) 



2'kR 

is the distance between tracers along the circumference of the ring in units of tracer size, 
and 

is a strongly temperature dependent prefactor. Fig.[5]shows that P is approximately constant 
for temperatures below 2 K but it increases sharply above 2 K. This means that, for the 



same geometry (ring's radius and number of tracers along the ring), the superfluid Stokes 
number is bigger for T > 2 K, hence the motion of vortex rings is much more disturbed by 
the presence of tracer particles than shown for example in Figs. [T]131 



IV. APPLICATION TO QUANTUM TURBULENCE 

Isotropic homogeneous quantum turbulence is characterized by the vortex line density 
L = A/V, where A is the total vortex length and V is the volume. In the first approximation 
this form of turbulence can be generated by applying a heat flux From the value of 
L one infers that the average vortex separation and the typical radius of curvature are of 
the order of £ ~ L^^/^. In a recent experiment [l^, solid hydrogen particles were used to 
visualize such turbulence: some tracers moved along the normal fluid in the direction of 
the heat flux, and other tracers were trapped on the quantized vortices as the vortex tangle 
drifted in the opposite direction. At T = 2 K the turbulence generated by the heat flux 



Q = 90 mW/cm^ corresponds to approximately 1^ £ ~ 0.006 cm. The hydrogen volume 
fraction used in the experiment was (p = Vhz/V ~ 10~^. Generalizing from a vortex ring of 
length 27iR to a vortex tangle of length A = LV and using = 4:N7ia^/3, the superfluid 
Stokes number can be estimated as 



30/5 fi 
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where 6 = Na/{2'nR) ^ {N/V){V/A)a = {N/V){a/L). If one tenth of the tracers are 
trapped inside vortices, using /3 ~ 20 (see Fig. [5]) we obtain 5* ~ 5 x 10^^ <^ 1. We conclude 
that in this particular experiment the tracers have not disturbed the dynamics of the vortex 
tangle. 



V. DISCUSSION 



We have derived the equations of motion of a vortex ring which contains buoyant tracer 
particles trapped on the vortex core. We have shown typical solutions of these equations to 
illustrate the difference between the motion of a bare, undisturbed vortex ring, and that of 
a vortex ring which is visualized by tracer particles. 

A flrst approximate measure of the disturbance caused by trapped particles is the geo- 
metrical quantity 6 given by Eq. (|27|) . the separation between tracers along the ring in unit 
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of tracer's size. Clearly 5 must be small for tracers not to disturb the vortex. A more pre- 
cise measure is the superfiuid Stokes number 5* = 5/5, where the dimensionless temperature 
dependent quantity /? is approximately constant for T < 2 K but becomes very large for 
T > 2 K. 

Finally, we have shown how this criterion is relevant to quantum turbulence, and can 
be used to determine the maximum hydrogen volume factor which can be used to visualize 
the turbulent flow without disturbing it. We have applied the criterion to a recent experi- 
ment [l^ in which the drift of vortex lines generated by a heat flux was measured and found 
that the disturbance was negligible. It would be interesting to apply the same criterion to 



experiments in which vortex reconnections were detected [l4|: this case is more challeng- 
ing, as during the reconnection process much smaller values of the radius of curvature are 
expected to arise which may be disturbed by tracer particles. 

VI. ACKNOWLEDGEMENTS 

We are grateful to W. F. Vinen. S. W. Van Sciver, M. S. Paoletti, D. P. Lathrop, and 
K. R. Sreenivasan for fruitful discussions. 



8 



[1] R. J. Donnelly, Quantized vortices in He //(Cambridge University Press, 1991). 
[2] P. M. Walmsley and A. I. Golov, Phys. Rev. Lett. 100, 245301 (2008). 

[3] P.-E. Roche, P. Diribarne, T. Didelot, O. Prangais, L. Rousseau, and H. Willaime, Europhys. 
Lett. 77, 66002 (2007). 

[4] D. I. Bradley, S. N. Fisher, A. M. Guenault, R. P. Haley, S. O'Sullivan, G. R. Pickett, and 

V. Tsepelin, Phys. Rev. Lett. 101, 065302 (2008). 
[5] T. V. Chagovets, A. V. Gordeev, and L. Skrbck, Phys. Rev. E 76, 027301 (2007). 
[6] V. B. Eltsov, A. I. Golov, R. de Graaf, R. Hanninen, M. Krusius, V. S. Lvov, and R. E. Sol- 

ntsev, Phys. Rev. Lett. 99, 265301 (2007). 
[7] W. F. Vinen and J. J. Niemela, J. Low Temp. Phys. 128, 167 (2002). 
[8] C. F. Barenghi, Physica D, 237 2195 (2008). 

[9] R. J. Donnelly, A. N. Karpetis, J. J. Niemela, K. R. Sreenivasan, W. F. Vinen, and 
C. M. White, J. Low Temp. Phys. 126, 327 (2002); D. Celik and S. W. Van Sciver, Exp. 
Therm. Fluid Sci. 26, 971 (2002); T. Zhang, D. Celik, and S. W. Van Sciver, J. Low Temp. 
Phys. 134, 985 (2004); T. Zhang and S. W. Van Sciver, J. Low Temp. Phys. 138, 865 (2005). 
T. Zhang and S. W. Van Sciver, Nat. Phys. 1, 36 (2005). G. P. Bewley, D. P. Lathrop, and 
K. R. Sreenivasan, Nature (London) 441, 588 (2006). 
[10] M. S. Paoletti, R. B. Fiorito, K. R. Sreenivasan, and D. P. Lathrop, J. Phys. Soc. Japan 77, 
111007 (2008). 

[11] M. S. Paoletti, M. E. Fisher, K. R. Sreenivasan, and D. P. Lathrop, Phys. Rev. Lett. 101, 
154501 (2008). 

[12] G. P. Bewley, M. S. Paoletti, K. R. Sreenivasan, and D. P. Lathrop, Proc. Nat. Acad. Science 
105, 13707, (2008). 

[13] S. Z. Alamri, A. J. Youd and C. F. Barenghi, Phys. Rev. Lett. 101, 215302 (2008). 

[14] D. R. Poole, C. F. Barenghi, Y. A. Sergeev, and W. F. Vinen, Phys. Rev. B 71, 064514 (2005). 

[15] Y. A. Sergeev, C. F. Barenghi, and D. Kivotides, Phys. Rev. B 74, 184506 (2006); erratum, 

ibid. 75, 019904(E) (2007). 
[16] D. Kivotides, C. F. Barenghi, and Y. A. Sergeev, Phys. Rev. B 75, 212502 (2007); ibid. 77, 

014527 (2008). 



9 



[17] C. F. Barenghi, R. J. Donnelly, and W. F. Vinen, J. Low Temp. Phys. 52, 189 (1983). 

[18] J. T. Tough, in Progress of Low Temperature Physics, vol. VIII, edited by D. F. Brewer (North 

Holland, 1982), pp. 133-219. 
[19] The heat flux Q induces the counterflow velocity Vns = Vn — Vs = Q/{PsS'T) where S' is the 

entropy. At T = 2 K S' = 0.956 J/(g K) and ps = 0.065 g/cm^, so Vns = 0.73 cm/s. The vortex 



line density is thus L = 7 ~ 2.7 x 10 cm , where, following Ref. 



13], 7 220 s/cm^ 



10 




I \ \ \ \ \ 1 

0.2 0.4 0.6 0.8 1 1.2 

t 

FIG. 1: (Color online) Dimensionless radius R of the vortex ring vs time t for: "a", solid (red) 
line: bare vortex ring (A^ = 0); "b", long-dashed (green) line: vortex ring with = 10 tracers, 
e = 0.19 X 10-^ C = 0.34 x lO^^; "c", short-dashed (blue) line: iV = 20, e = 0.37 x IQ-^, 
C = 0.68 X 10-^ "d", dotted (magenta) line: iV = 40, e = 0.75 x 10-^ C = 0.14; "e", dashed- 
dotted (cyan) line: AT = 80, e = 0.15 x 10-^ C = 0.27. 
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FIG. 2: (Color online) Dimensionless radial velocity u of vortex ring vs time t for = (bare 
ring, solid line labeled "a") and = 10, 20, 40 and 80 tracers. (Labels and line styles (colors) 
correspond to those of Fig. [TJ) 
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FIG. 3: (Color online) Dimensionless axial velocity v of vortex ring vs time t for = (bare ring, 
solid (red) line labeled "a") and N = 10, 20, 40 and 80 tracers. Labels and line styles (colors) 
correspond to those of Figs. [1] and [2j 
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FIG. 4: (Color online) Superfluid Stokes number vs dimensionless time t for N = 10, 20, 40 and 
80 tracers. Labels and line styles (colors) correspond to those of Figs. [l][3l 
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